Intrinsic finite-size effects in the two-dimensional XY model 
with irrational frustration 
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This study investigates in detail the finite-size scahng of the two-dimensional irrationally frus- 
trated XY model. By means of Monte Carlo simulations with entropic sampling, we exaiuine the 
size dependence of the specific heat, and find remarkable deviation from the conventional finite-size 
scaling theory, which reveals novel intrinsic finite-size effects. Relaxation dynamics of the system 
is also considered, and correspondingly, finite-size scaling of the relaxation time is examined, again 
giving evidence for the intrinsic finite-size effects and suggesting a zero-temperature glass transition. 
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The two-dimensional XY model with uniform frustra- 
tion has attracted much interest in connection with quite 
a few physical realizations such as periodic arrays of cou- 
pled Josephson junctions in uniform magnetic fields [y . 
Of particular interest is the case of irrational frustration, 
where the possibility of a structural glass phase with- 
out intrinsic random disorder has been a long-standing 
question. In the model with irrational frustration, the 
periodicity of the ground state should be incommensu- 
rate with the underlying array periodicity and the ground 
states are infinitely degenerate, suggestive of the absence 
of a finite-temperature transition. It was thus argued 
that there is no phase transition at finite temperatures 
in the thermodynamic limit H , which was also favored by 
numerical simulations of the current-voltage {IV) char- 
acteristics pi. The latter, performed with resistively 
shunted junction (RSJ) dynamics, demonstrated that the 
IV scaling is consistent with a zero-temperature vortex- 
glass transition proposed in Ref. B. In other simulations, 
on the other hand, a metastable finite-temperature glass 
transition was proposed [^jj^] and the similarity of the 
relaxation dynamics to the primary relaxation in a su- 
percooled liquid was pointed out 0. Subsequent Monte 
Carlo (MC) simulations of the original XY model [g| and 
of the corresponding Coulomb gas 0] then suggested the 
relevance of boundary conditions and dynamics in ob- 
serving the reported glass effects, leaving the existence 
of a finite-temperature transition still inconclusive. 

One way to resolve this is to attribute the observation 
of the finite transition temperature to intrinsic finite- 
size effects. In a finite system the effects of irrational 
frustration may not be distinguished from those of ratio- 
nal frustration, the value of which is given by a ratio- 
nal approximant sufficiently close to the irrational pO[ . 
Consequently, unlike the system in the thermodynamic 
limit, a finite system (with irrational frustration) tends 
to behave as if the value of the frustration were given by 
the rational approximant. Accordingly, the finite system 



may display such transition-like behavior as the specific 
heat peak (rounded off due to usual finite-size effects) 
at a finite temperature, even though the system in the 
thermodynamic limit undergoes a zero-temperature tran- 
sition, as argued in Ref. [0 . Here the approximant which 
properly describes a finite system depends on the sys- 
tem size [nol , raising the interesting possibility of intrin- 
sic finite-size effects in addition to the usual ones; this 
has not been probed in previous simulations, where, for 
simplicity, fixed rational approximants were mostly used 
instead of the actual irrational value of frustration. 

In this Letter, we adopt the idea of the intrinsic finite- 
size effects to resolve the controversy on the phase tran- 
sition in the two-dimensional irrationally frustrated XY 
{\FXY) model. For this purpose, we use the exact irra- 
tional number up to the machine precision, and perform 
extensive MC simulations with improved entropic sam- 
pling |ll| to compute the specific heat. The obtained 
data turn out to be inconsistent with the conventional 
scaling theory of either the first-order transition or the 
second-order one, revealing the presence of novel intrin- 
sic finite-size effects. To confirm the existence of such in- 
trinsic finite-size effects, we also study the relaxation dy- 
namics of the system and construct the finite-size scaling 
theory of the relaxation time, on the basis of the zero- 
temperature glass transition together with the intrinsic 
finite-size effects. The size dependence of the relaxation 
time indeed agrees well with the scaling theory, providing 
strong evidence for the intrinsic finite size effects. 

The Hamiltonian for the irrationally frustrated XY 
model on an L x L square lattice is given by 



H 



-jyj cos( 



A 



(1) 



where the angle 0i corresponds to the phase of the super- 
conducting order parameter at site i and J represents the 
Josephson coupling strength in the case of a Josephson- 
j unction array. The plaquette sum of the bond angle Aij 



is constant over the whole lattice, ^ Aij — 27r/, where 
the frustration parameter / is set equal to the golden 
mean: f ^ fl = {V5 - l)/2. 

To investigate the finite-size effects of this system, we 
consider systems of size L — 5,8, 13, 2f, and 34, the en- 
tropies of which are computed by means of extensive 
MC simulations with entropic sampling. The entropy 
obtained for each system allows one to calculate the spe- 
cific heat of the system. Both periodic boundary con- 
ditions (PBC) and free boundary conditions (FBC) are 
used and the results are compared with each other, from 
which two types of finite-size effects are identified. In 
principle both boundary conditions should yield correct 
limiting behavior as the system size L is increased, al- 
though for finite L the PBC induces mismatch at the 
boundary. The values of L are thus chosen in such a way 
that the mismatch present in the PBC is kept minimum. 
In the case of the golden mean, it can be achieved by 
choosing the Fibonacci numbers. 

Figure presents the temperature dependence of the 
obtained specific heat in the PBC for various system 
sizes. For L = 34, in order to get reliable results, we 
have performed as many as 7.5 x 10^° MC steps per site 
using a Cray T3E supercomputer system. As L gets in- 
creased, the peak is shown to grow and the transition 
region to narrow whereas the peak position shifts toward 
lower temperatures. The position of the specific heat 
peak in the system of size L, denoted by Tm{L), is ex- 
pected to approach the critical temperature Tc as L is 
increased. The well-known finite-size scaling theory of 
the second-order transition [ p2|p3| gives the behavior 

T^{L)-T,^L-^'\ (2) 

to which the least-square fit of our data leads to Tc = 
and 1/v = 0.36. On the other hand, the peak height 
in Fig. |l|, expected to grow as L°''^, gives the value 
a/v w 0.17, which, via the Josephson scaling relation 
{2 ~ a — dv with the system dimension d = 2), leads 
to the value v k, 0.92; this is consistent with the value 
u « 0.9 obtained in the numerical study of IV charac- 
teristics g. Thus the result of the fit manifests inap- 
plicability of the conventional finite-size scaling theory. 
Further, it is also inconsistent with the scaling theory of 
the first-order transition ||l^,n, according to which 1/v 
in Eq. (0) is replaced by the system dimension d — 2. It is 
thus concluded that the system does not display conven- 
tional finite-size effects associated with either a first-order 
transition or a second-order one. Still remains the possi- 
bility of a weakly first-order transition, in which case the 
system may appear to exhibit size dependence similar to 
that for a second-order transition |l3||l^ ]. Such pseudo- 
divergent scaling behavior should disappear as the system 
size exceeds the correlation length of the system. 

Deviation from these conventional (second-order, first- 
order, or weakly first-order) transitions can be observed 
most clearly in comparing the size dependences of the 
specific heat under PBC and under FBC. Periodic bound- 
ary conditions tend to suppress fluctuations, exhibiting 



the peak position at a temperature higher than the true 
critical temperature T^. Under FBC, on the other hand, 
no such suppression exists and large fluctuations in gen- 
eral yield the peak at a lower temperature [|l2[. Thus 
the ordinary finite-size effects are manifested by the shift 
of the peak position toward Tc as the system size is 
increased, from higher temperatures in PBC and from 
lower temperatures in FBC, respectively. The results of 
test runs for the fully frustrated XY model with / — 1/2, 
shown in Fig. 0(a), indeed demonstrate such behavior of 
the ordinary finite-size effects: The data in both cases, 
PBC and FBC, follow Eq. (|) (up to the overall sign), and 
the least-square fit with the values of Tc and v in Rcf . Ill] 
yields the lines displaying nice scaling behaviors ||lj . On 
the other hand. Fig. ||(b) reveals quite different features 
of the finite-size effects in the IFXY model. Here the 
peak position Tm{L) in PBC is sensitive to the system 
size LjTCducing rapidly with L, while it is not the case in 
FBC iQ. In particular, unlike the case of rational frus- 
tration, the peak position in FBC does not shift toward 
higher temperatures, violating the conventional scaling 
given by Eq. (||). 

Such peculiar behavior can be explained in terms of 
intrinsic finite-size effects present in addition to the or- 
dinary ones pQ l : For the frustration given by the golden 
mean, the appropriate rational approximants are given 
by the ratios between adjacent pairs in the Fibonacci se- 
quence, i.e., fk = qk~i/qk, where qk denotes the fcth Fi- 
bonacci number. Among these rational frustration values 
/fc's, those satisfying the following condition determine 
nature of the transition in the (finite) system of size L: 



<?fc<L< \f-fkl 



(3) 



where f — il. While the first inequality in Eq. (y) sim- 
ply implies that the system size should be sufficiently 
larger than the size of a unit cell, the second one, related 
with the charge neutrality condition, gives the criterion 
that the defects due to the mismatch f — fk are negligible 
since the corresponding phase change across the system is 
27rL(/ — fk) |J0[. Accordingly, as the system size L is in- 
creased, the mismatch should get smaller and the bigger 
Fibonacci number qk is needed to describe the transition 
at the size. The transition temperature as well as nature 
of the transition depends on the frustration parameter or 
on qk ||3,0,ll^ and thus on the system size. 

It is therefore plausible to replace Tc in Eq. (0) by 
Tc{L), the transition temperature of the system of size 
L, which yields 



T„(i) - Tc{L) = Cj^L-^'" 



T,n{L) - Tc[L) = -Cf 



^-1/. 



(4) 



with appropriate (positive) amplitudes Cp and Cf for PBC 
and for FBC, respectively. Here Tc{L) is expected to de- 
crease with the system size, approaching the true critical 



value Tc in the thermodynamic hmit iQ • Since the pre- 
cise hmiting behavior is unknown, we assume the simple 
behavior 



Tc{L) -Tc ^ CiL 



(5) 



Note that in contrast to the ordinary finite-size effects, 
which depend crucially upon the boundary conditions, 
the intrinsic effects described by Eq. (g) change the tran- 
sition temperature itself, regardless of the boundary con- 
ditions. 

Equations (0) and d|) thus yield the appropriate scal- 
ing form for the IFXF model: In PBC, both ordinary 
and intrinsic effects add up to give the finite-size scaling 
in the form 



T„,{L) - T, - qL-" + CpL-l/^ 



(6) 



whereas in the case of FBC, the two effects give the con- 
vergence behavior in the opposite direction, leading to 
the form 



T^{L)-Tc^c,L-"^-CfL-^'\ 



(7) 



Indeed using the value v = 0.9, we can fit the peak po- 
sition data both in PBC and in FBC io Eqs. (P) and ffl) 
with Tc = 0, respectively, to obtain the same value 



0.28 ±0.05. 



(8) 



The resulting scaling curves are displayed in Fig. g(b), 
giving strong evidence for Tc = 0. 

To confirm this, we also study the relaxation dynamics 
of an L X L RSJ array under PBC in both directions p9[ . 
To obtain dynamical behavior, we integrate directly the 
equations of motion with the time step Ai = 0.05, start- 
ing from random initial configurations and taking aver- 
ages over 300 to 30000 independent runs. In the limit 
i — !■ c», the energy function is expected to approach the 
equilibrium value at given temperature T, allowing us 
to estimate the relaxation time r of the energy function. 
The obtained behavior of r is shown in Fig. for various 
sizes L = 5,8, 13,21, and 34. Notice that for T > 0.25 
the data do not depend appreciably on the system size, 
which explains why the size effects were not observed in 
Ref. g. Still the scaling of IV characteristics, not rely- 
ing much on the system size in the temperature region 
probed, can give an accurate value of the exponent v [^. 

To examine the finite-size effects in the relaxation time 
data, we construct the suitable scaling form of the re- 
laxation time. In the presumed zero-temperature glass 
transition, the relaxation has an activated form, and the 
relaxation time diverges exponentially as the tempera- 
ture approaches zero, taking the form ll,|2C 



T = Aexp(A£;/T). 



(9) 



Here AE is the typical barrier that a vortex should over- 
come to move across the correlation length ^, and scales 



as AE ^ ^'^ with the barrier exponent ip. In the pres- 
ence of intrinsic finite-size effects, which lead to the size- 
dependent transition temperature Tc{L), the correlation 
length should display the scaling behavior 



C = T/(iV-[T-Tc(T)]) 



(10) 



with the appropriate scaling function f{x) such that 
f{x — > 0) — > constant and f{x -^ oo) -^ x~'^ . The cor- 
responding finite-size scaling of the relaxation time thus 
takes the form 



logT = i^(Li/''[T-Tc(T)]) 



(11) 



with the scaling function g{x) ex [f{x)]'^. 

Figure ^ presents the scaling plot of the relaxation 
time T of the IFXY model, showing that the data in- 
deed fit well to the finite size-scaling form in Eq. ([ll|). 
The critical temperature Tc{L) at each size L has been 
estimated from the specific heat data. /^From this scaling 
with -0 = 0.54, we obtain the exponent v — 0.89 ± 0.2, 
which is in good agreement with the results of the specific 
heat data and of the IV characteristics in Ref. [^ . 

In summary, we have studied the finite-size scal- 
ing of the specific heat and the relaxation time in 
the two-dimensional irrationally frustrated XY model. 
This has revealed intrinsic finite-size effects and a zero- 
temperature glass transition in the thermodynamic limit, 
thus resolving the controversy as to the existence of a 
finite-temperature transition. Since the analytic argu- 
ment for the intrinsic effects is rather general pX| ], we 
believe the conclusion to be valid for any irrationals; this 
is indeed supported by the preliminary results for the sil- 
ver mean case, the detailed investigation of which is left 
for further study. Finally, we point out that for generic 
values of frustration, a quasiperiodic array displays be- 
havior similar to that of a periodic array with irrational 
frustration |2^ . This suggests the presence of intrinsic 
finite-size effects also in the quasiperiodic array, as mani- 
fested by the apparent glass behavior reported in Ref. [|| . 
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FIG. 1. Specific heat versus temperature in the irrationally 
frustrated XY model under periodic boundary conditions, for 
size L = 5, 8, 13, 21, and 34. The peak position shifts toward 
lower temperatures as the system size is increased. 
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FIG. 2. Size dependence of the position of the specific heat 
peak, (a) In the rational case (/ = 1/2), the data in both 
PBC and FBC follow Eq. (S, with the lines denoting the 
least square fits to Eq. (^. (b) In the irrational case (/ = f2), 
the data fit to Eq. (|) (PBC) or (0) (FBC) rather than to 
Eq. (m). In (a) error bars are not larger than the symbols. 
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FIG. 3. Relaxation time r versus the temperature T in the 
semi-log scale for system size L = 5, 8, 13, 21, and 34. Error 
bars have been estimated from the standard deviation and 
those data points without them have errors not larger than 
the size of the symbols. The dotted lines are merely guides 
to the eye. 
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FIG. 4. Fit of the relaxation time data to the fi- 
nite-size-scaling formula given by Eq. ([l_l|), with ip = 0.54 
and 1/u — 1.12. The critical temperature Tc(L) has been 
obtained from the specific heat data. 



